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^ ' Abstract. We present a new approach to study measures on ensembles of contours, polymers 
O 

or other objects interacting by some sort of exclusion condition. For concreteness we develop 

O . 

OO ' it here for the case of Peierls contours. Unhke existing methods, which are based on cluster- 



expansion formalisms and/or complex analysis, our method is strictly probabihstic and hence 
can be applied even in the absence of analyticity properties. It involves a Harris graphical 
construction of a loss network for which the measure of interest is invariant. The existence of the 



^ process and its mixing properties depend on the absence of infinite clusters for a dual oriented 
Ci \ percolation process which we dominate by a multitype branching process. Within the region of 
subcriticality of this branching process the approach yields: (i) exponential convergence to the 
equilibrium measure, (ii) clustering and finite-effect properties of the contour measure, (iii) a 
particularly strong form of the central limit theorem, and (iv) a Poisson approximation for the 
distribution of contours at low temperature. 
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1 Introduction 

Contours were introduced by Peierls (1936) to prove the existence of a first-order phase tran- 
sition for the Ising model in 2 or more dimensions. His argument, later put on a rigorous 
mathematical basis by Dobrushin (1965) and Griffiths (1964), used contours only as an auxil- 
iary device to estimate spin correlations. Polymer models, in the sense of interest here, were 
introduced later by Gruber and Kunz (1971). These are abstract general models of possibly 
extended objects that interact only by volume exclusion. They include both contour ensembles 
and ensembles formed by the open walks ("polymers") or surfaces obtained in high-temperature 
expansions. Gruber and Kunz (1971) were the first to treat these models as probability ensem- 
bles of their own, and to ask genuinely probabilistic questions such as existence and properties 
of the corresponding probability measure. The formalism of cluster expansions, whose use in 
mathematical physics started with a paper by Glimm, Jaffe and Spencer (1976), soon estab- 
lished itself as the technique of choice to study these type of systems and questions [Malyshev 
(1980), Seller (1982), Brydges (1984)]. The formalism was extended by Kotecky and Preiss 
(1986) to objects obeying generalized exclusion laws defined by compatibility relations. This 
extension was taken up by Dobrushin (1996, 1996a) who proposed to call animal models to 
such general systems and introduced a new approach to the construction of the expansions. 

The cluster-expansion technology, being designed to construct and study distributions of 
general systems with exclusions, seems potentially very useful for probabilists in general. Nev- 
ertheless, its use has so far remained confined to the mathematical physics community work- 
ing in statistical mechanics and quantum field theory. This unfortunate situation has been 
pointed out by Dobrushin (1996a), who attributed it to two reasons: (1) "its analytical and 
combinatorial complexity" , and (2) "the absence ... [of a] systematical exposition oriented to 
mathematicians". He addressed both issues in his posthumous review, Dobrushin (1996a), 
where he presented an exposition geared towards "probabilistic interpretations and applica- 
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tions" , based on his new approach that avoids "tremendous combinatorial considerations" . In 
fact, his approach does not resort to cluster expansions at all. 

In our opinion there is still an additional aspect that explains the lack of popularity, among 
probabilists, of this powerful technique: All the existing formulations transcend the probabilistic 
framework. First, the expansions used are a bit unnatural from the measure-theoretical point 
of view. Cluster expansions were, in fact, originally introduced to control the pressure of 
gases with exclusions. The rigorous proof that they converge and have nice mathematical 
properties required highly nontrivial combinatorial estimates which took a reasonable form 
only with the insight of Cammarotta (1982) [clearly described in Brydges (1984) and Pfister 
(1991)]. The existence of a measure is proven by using pressure-like expansions for numerator 
and denominator and cancelling out terms. It is, therefore, a rather indirect approach whose 
mathematical bottleneck refers to an object — the pressure — that from the probabilistic point 
of view is just auxiliary. Dobrushin's approach, on the other hand, avoids the use of the 
pressure and the cluster expansion but its use of complex analysis reveals that the corresponding 
hypotheses and results go beyond probability. Furthermore, for actual computations one needs 
to go back to the traditional approach and its explicit expressions for the correlation functions 
(or the pressure). 

A second manifestation of "probabilistic unnaturalness" is conveyed by the results them- 
selves. Indeed, the existing formulations require the absolute convergence of the expansions 
involved. As a consequence, besides existence and mixing properties, they yield analyticity of 
the correlation functions with respect to different parameters, for instance with respect to the 
exponential of minus the inverse temperature. Though analyticity is a very nice property to 
have — in particular it allows Dobrushin to produce an amazingly simple proof of the central 
limit theorem — it is also a symptom that these approaches are too strong and not optimal 
from the probabilist point of view. This is not just an academic remark. The most interesting 
recent applications of cluster-expansion methods fall outside these formulations, as they involve 
measures that are known or suspected to have non-analytical behavior: Measures at intermedi- 
ate temperatures [Olivieri (1988), Olivieri and Picco(1990)], measures for annealed disordered 
systems [von Dreifus, Klein and Perez (1995)], measures for long-range interactions [Bricmont 



3 



and Kupiainen (1996)] and infinite-dimensional Sinai-Ruelle-Bowen measures [Bricmont and 
Kupiainen (1997)]. 

In this paper we present a novel approach to the study of animal models which presents 
a number of advantages regarding these issues. For concreteness we discuss the case of usual 
Peierls contours; a more general treatment will be presented in Fernandez, Ferrari and Gar- 
cia (1998a). Here is an overview of the main features of our approach. 

1. The approach is purely probabilistic, no cluster expansion or complex analysis is involved. 
The measure is obtained as the unique stationary measure of a Markov process. The 
condition of validity of our theory is stated in terms of a backwards oriented percolation 
process. The theory holds when percolation is absent. 

2. The range of validity of the theory exceeds that of previous approaches [see comment 
after ( |2.10D ]. Within this range we obtain all the properties yielded by the latter — 
existence, uniqueness, exponential mixing, central limit theorem — with one conspicuous, 
and expected, exception: analyticity. 

3. We obtain a rather nice version of the central limit theorem [stronger than that in Do- 
brushin (1996a)]. 

4. The approach allows us to show that the rescaled distribution of contours of a fixed 
length convergence towards a Poisson process. We are not aware of similar results in the 
literature. 

5. The construction constitutes, in fact, a simulation scheme that converges to equilibrium 
exponentially fast. Hence, it has the potential to become a very efficient computational 
tool. 

6. The avoidance of series expansions for the pressure makes our approach more direct to 
compute general properties of the equilibrium measure, but limits its use for the esti- 
mation of "thermodynamic" quantities. For instance, the approach does not seem to be 
suitable for the study of "surface corrections" to the presure. Bounds on these corrections 
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arc crucial for several applications of contour ensembles [see eg. Zahradnick (1984), Borgs 
and Imbrie (1989)]. 

In this paper we present a careful statement of these results and a sketch of their proofs. 
We aim at providing a streamlined exposition free of inessential technicalities that may obscure 
the natural form of the construction. Nevertheless, we present enough details for an educated 
probabilist to reconstruct most of the missing links. The full argument will be presented in 
Fernandez, Ferrari and Garcia (1998), theretofore referred as FFG. 

2 Contour distribution and loss networks. Results. 
2.1 Contours 

The contours for the ferromagnetic Ising model with "+1" boundary conditions, in dimensions 
d > 2, are surfaces constructed with {d — l)-dimensional unit cubes — traditionally known 
as plaquettes — centered at points of Z'^ and perpendicular to the edges of the dual lattice 
Z'^ + (i, • • • , |). We shall identify a plaquette with its center and denote x e 7 if the plaquette 
centered at x is in 7. Two plaquettes are adjacent if they have a common [d — 2)-dimensional 
face. A collection of plaquettes forms a connected surface if for every two plaquettes x, y one 
can find a finite sequence of plaquettes, starting at x and ending at y, such that two consecutive 
plaquettes of the sequence are adjacent. A closed surface has every {d — 2) -dimensional face 
shared by 2 or 4 plaquettes. A contour, 7, is a connected and closed family of plaquettes. We 
say that two contours 7 and 9 are incompatible, and denote 7 fl ^ 7^ 0, if they have adjacent 
plaquettes. We use the notation \x — y\ for the minimal number of plaquettes needed to link, 
in a connected fashion, x with y (this is also known as "Manhattan distance"). 

For A C 7/, denote by ^(A) the set of contours whose plaquettes have centers in A. A 
configuration of contours r] G N^(^) is a function that at each contour 7 assigns a natural 
number 77(7) indicating the number of contours 7 present in tj. The subset A" (A) C N^(^) of 
compatible-contour configurations is defined as 

X{A) = {r]e {0, 1}^(^) ; 7^(7) r]{e) = if 7 n ^ 7^ 0} (2.1) 
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that is, a configuration of contours is compatible if it contains at most one copy of each contour 
and does not contain two intersecting contours. 

For each fixed P e M+, a parameter usually called the inverse temperature and for each 
finite A C Z*^ define the measure //^ on X{A) by 

Av) = ^ ^ ^ (2.2) 

where I7I is the area (^number of plaquettes) of the contour 7 and is a normalization 
constant making a probability. 

2.2 Loss network of contours 

We introduce a birth-and-death dynamics on the set of compatible contours. This process is 
known in the literature as loss network, see Kelly (1991) and references therein. 
We define the process rj^ as a Markov process on A" (A) with generator given by: 

A^m= E e-^^Mv^^ ^ ximnv^') - m]+ E ^(7)[/(r")-/(^)] (2.3) 

■ye6(A) 7Ge(A) 
for / : X{A) — > R, where 1{ • } denotes the characteristic function of the set { • } and for 

( 7](e) if ^ 7^ 7 

[ rj{-i) ±1 if ^ = 7 
It is immediate to check that the measure /x^ is reversible for r]^. 

In terms of loss network language, the above process can be described as follows. Consider a 
network consisting of a finite number of links represented by plaquettes with vertices in A C Z*^, 
each link comprising one circuit. Calls are offered to this network along routes 7 G Q{A) 
according to independent Poisson processes with rate e^^'''''. A call accepted on route 7 holds 
all links along this route for an exponential holding time with mean 1 and on completion of 
the service releases all these circuits simultaneously. All arrival streams and holding times are 
mutually independent. A call is accepted along route 7 e ^(A) if 7 is not compatible with 
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other calls already in progress. Hence rj^ = {vtil))i<^G{A-) where rj^i'y) is the number of calls in 
progress on route 7 at time t, then a call is accepted along route 7 at time t if 

We can represent this model as a solution of the following system of equations: 

vHl)=V^il)+ fi{ E vtil') = 0}dN^ie-^^'^s)-N;[[\^i^)ds) (2.5) 

where iV+ e are independent unit Poisson processes; A^^ creates new contours and 
destroys them. 



2.3 Range of validity of the approach 



Let X = {r] e {0, l}^^^ ^ : vilM^) = if 7n^ ^ 0}. Since GiZ"^) is countable, X is compact in 
the product topology. Let / be a continuous function on X. The infinite- volume loss network 
on X has formal generator given by 

Afiv) = Y.e-f'^Mv^^ e X}[f{r^^r) - /(r^)] + ^r^(7)[/(r^-^) - /(r^)] (2.6) 

y&G 



where rj^'' was defined in (p.4|) 



We use a graphical construction to show that a sufficient condition for the existence of a 
process r]t on X with generator A is 

= ^ I7I e-^l^l < 00. (2.7) 
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Using the fact that X is compact, abstract nonsense imply that, under ( |2.7|) , there exists an 
invariant measure fi for rit. However, the way of proving existence is so general that we are 
not able to show any further property of this measure. We remark that, as pointed out by 
Aizenman, Bricmont and Lebowitz (1987), ( |2.7| ) defines a "Peierls" inverse temperature, 

/5p = inf{/3 : < 00} , (2.8) 



above which, with probabihty one, only a finite number of contours surround any given site (a 
fact that, for the Ising model, implies existence of spontaneous magnetization). The results of 
this paper, however, apply to the more limited regime 

13 > Pm, (2.9) 

where 

/3M = inf{/5:A^<l/(d-l)}. (2.10) 

For the Peierls contours the best estimations of the range of validity of "traditional" cluster- 
expansion approaches follow from Proposition 5.6 in Dobrushin (1996a), which has been stated 
in its most precise form by Lebowitz and Mazel (1997). As a matter of fact, these authors 
present their estimations in a form slightly different to ours: They consider contours with a 
given site of the dual lattice in its interior, rather than contours containing a given plaquette as 
we do. The final expressions obtained in these two cases are not directly comparable because 
they involve differently-aimed upper bounds. For a meaningful comparison we have either 
to transcribe our approach in terms of interior sites, or to write theirs in terms of anchoring 
plaquettes. The latter policy leads to a bound 



^e^l^le-^l-'l < (2.11) 



d - 

730 

for some constant c > 0. This bound can be read off the work of Lebowitz and Mazel (1997) 
[who obtain c = /3e~°'^''^], where in fact all the hard estimates [from their formula (2.7) till 
the end of their paper] refer to contours containing a fixed plaquette. On the other hand, our 
condition (|2.9|) - (|2.1CI|) implies 

El^l^"""'^ rf3T' (2.12) 

790 

which is strictly weaker than (|2.11|) because > x for x > 1. 

Lebowitz and Mazel show that, defining (3lm as the infimum of (3 satisfying ( |2.11| ), 

(3lm > 64^^, (2.13) 

where ( |2.12| ) plus their counting method, yields 

/3m > 6^. (2.14) 
a 
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On the other hand, Aizenmann, Bricmont and Lebowitz (1987) show that the Peierls temper- 
ature defined by ( |2.8| ) satisfies 

,.>J% (2.5, 

These three temperatures mark, therefore, hmits where different properties can be proven by 
perturbation arguments. For (3 > Pp, each site of the dual lattice is sorrounded by a finite 
number of contours. In spin language, this means lack of percolation of minority spins (which, 
in turn, implies symmetry breaking and, by FKG, non zero magnetization). For f3 > /3m, in 
addition, properties Rl — R5 listed below can be proven by cluster-expansion-like methods. 
Finally when /? > Plm methods of this type also yield analytic temperature dependence. 



2.4 Results 

We say that / has support in T C Z'^ if / depends only on contours intersecting T (not 
necessarily contained in T). Let |Supp(/)| = min{|T| : / has support in T}. When we write 
Supp (/) we mean any T such that |T| = |Supp (/)| and / has support in T. For instance, if 
fiv) = ^7(7)) Supp (/) may be set as {x} for any x G 7. 

A closer analysis of the graphical construction allows us to show that for /5 > Pm the 
following results hold. These are our main results. 

Rl. Reversibility and uniqueness: there exists a unique invariant measure n for rjt. Further- 
more, II is reversible for the process rjt. 

R2. The rate of convergence to the invariant measure is exponential. Let S^S{t) be the distri- 
bution of the process at time t when the initial configuration is ^. For measurable / we 
prove 

sup\i,f-S^S{t)f\ < ||/||oo|Supp(/)|e-^°* (2.16) 
for any Mq < {1 - (d - l)A^)/(2 - {d - 1)A^). 

R3. Infinite- volume limit: Let A be a (finite or infinite) subset of Z"^ and / a measurable 
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function depending on contours contained in A. Then 

|/i/-/iVl<ll/llooM2 Yl e-^3<.,A^) (2.17) 

3::gSupp(/) 

where = e^^-^^)/^'^-^) and M3 > {P-pM)/{d-l). We denoted d{x, A'^) = mm{\x-y\ : 

y e A^}. 

R4. Clustering. For measurable functions / and g depending on contours contained in an 
arbitrary set A C Z"^: 



<2||/||oo|kl|oo(M2f \x-y\e-''-^\^~y\ (2.18) 

xGSupp(/), 
j/GSupp (g) 

where M2 ad M3 are the same of (|2.17|). This includes the infinite- volume measure 



R5. Central limit theorem. Let / be a measurable function on X with finite support such 
that fif = and /i(|/p+^) < 00 for some 6 > 0. Assume D = J2x ^^if'^^f) ^ ^- Then 
D < 00 and 

-^y^rj =^ NormalfO,^) (2.19) 

where the double arrow means convergence in distribution. This result generalizes (the 
central limit) Theorem 7.4 of Dobrushin (1996a). In the latter, only functions depending 
on a finite number of contours are considered. 

For the following result we write fip to stress the /? dependence of fi. 

R6. Poisson approximation. Let r]'^ distributed with /x^. For each measurable y C M'^ let 

V{a) = {xeZ'^ : x/a G V^}. (2.20) 

For each j, the process N^'^ defined by 

m^{v)= Yl ^''W- (2-21) 

7Cy(e«),|7|=i 

converges weakly to a unit Poisson process on M'^ as /3 ^ 00. The rate of convergence is 
exponential in /3. 
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The key to the proof of the above results is a graphical construction of the process starting 
from a marked Poisson process in Z'^ x M. The marks determine random cylinders whose bases 
are the contours and the heights are exponentially distributed random times. The exclusion 
condition is imposed through the study of the "ancestors" of each cylinder (Section These 
ancestors determine a (backwards) oriented percolation process, and our construction is feasible 
if there is no such percolation. This is the meaning of the condition (3 > I3m- All our results 
follow from the estimation of the spatial and temporal extension of the cluster of a (finite 
number of) cylinders(s) (Section ^ and This estimation is done through a domination of 
the percolation process by a multiple branching process which, in the regime (3 > (3m has 
exponential moments (Section ^). The proof of R6, also based in the above properties, is 
omitted here. A complete proof is presented in FFG. 

Ferrari and Garcia (1998) used space-time percolation to show ergodicity of loss networks 
under low arrival-rate of calls. 

3 Graphical construction. The BO-cluster 
3.1 Finite volume 

To each contour 7 G ^ we associate a Poisson process of rate e"^''''', and to each time event 
Tfc(7) of the Poisson process we associate an independent exponentially distributed time S'fc(7) 
of mean one. The collection C = (Tfc(7), 5'fc(7))^gg,fcgz is a family of double-sided independent 
marked Poisson processes, with the convention T__i(7) < < To (7). The /cth attempt of birth of 
a contour 7 occurs at time Tk{'~f)] 5*^(7) corresponds to the lifetime of the contour. Each triplet 
(7, Tfc(7), 5*^(7)) is called a cylinder of basis 7 birth-time Tfc(7) and lifetime Ski'y). To each 
contour 6 present in the initial configuration rjo = t] we independently associate an exponential 
time S{9) and cylinder (6*, 0, S{6)). The collection of initial cylinders is called C(0). We realize 
the dynamics (deterministic) function of C and C(0). 

When the number of possible contours is finite, the construction for t > is as follows. We 
construct inductively K[o,t], the set of kept cylinders. The complementary set corresponds to 
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erased cylinders. First include all cylinders of C(0) in Kp,*]. Then, move forward in time and 
consider the first Poisson mark: The corresponding cylinder is erased if it intersects any of the 
cylinders already in Kjo^tj, otherwise it is kept. This procedure is successively performed mark 
by mark until all cylinders born before t are considered. Define rjt G Af(A) as 

Vt{l) = r7o(7)l{5(7) >t} + l{3A; : (Tfc(7), Tfc(7) + ^^(7)) 9 t and (7, Tfc(7), 5^(7)) is kept}, 

(3.1) 

that is, rjt signals all contours which are basis of a kept cylinder that is alive at time t. 

It is tedious but easy to show that rjt has as generator an operator defined as A, but with 
the sums restricted to the finite set of contours involved. In particular, when the contours are 
contained in a finite region A, we obtain the process rj^ with generator A^. 

The above finite- volume construction can also be performed in (— cxd, 00). Indeed, rjf is an 
irreducible Markov process in a finite state space. Hence, with probability one there exists a 
sequence of ordered random times ^^(C) such that no cylinder in C is alive by time tfc(C). 
Furthermore IE(tfc_(_i — tfc) < 00. Therefore one can apply the above construction independently 
in each interval [tk{C),tk+i{C)). In this case the cylinders of C(0) play no role. This procedure 
is time-translation invariant and so is the distribution of rjt. This distribution is precisely given 
by the measure /i^. 

3.2 Infinite volume 

For infinite volume, the Poisson processes are indexed by an infinite set of contours. Hence, it 
is not possible to decide which is the first mark in time. The construction must be performed 
more carefully. There are two alternatives. 

The first alternative is to divide the time interval [0, t\ in successive intervals of small length 
h and perform the construction in each one of those intervals. Under ( p.7|) and for small h, it 
is possible to partition 7/" in finite regions such that each contour born in [0, h] is contained 
in exactly one of these regions. To show this, one considers the percolation of (projected) 
contours and dominates the area occupied by the contours by a branching process. Such a 
construction is at the heart of Harris (1972) original graphical construction of particle systems 
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and it is reviewed by Durrett (1995). The mark-by-mark construction described above can be 
performed in each of these finite regions to construct the process in the time interval [0, h]. The 
same procedure can then be apphed in the interval [h, 2h], etc. 

The second alternative is the one we really use. In order to know whether a cylinder C G C 
is kept, one has to look at the set of cylinders C (born before C and) alive at the birth-time 
of C whose basis intersects the basis of C. This set is called the first generation of ancestors 
of C . The second generation of ancestors of C consists, previsibly, of the ancestors of the 
ancestors, that is those cylinders that are in the first generation of ancestors of some C in the 
first generation of ancestors of C . Recursively we construct in this way the nth generation of 
ancestors of C . The set of ancestors (of any generation) of C is called the BO- (backwards 
oriented) cluster of C and it is denoted by A(C). This set may contain cylinders in C(0). [We 
remark that this BO-cluster is a cluster of space-time cylinders, it is different from the usual 
cluster of contours considered in the classical works on cluster expansions.] If for some C the 
BO-cluster of C has a finite number of cylinders, then we can decide whether C is kept or not 
by looking at A[o,i](C) = {C G A(C) : C is born in [0,t]}, the set cylinders in the BO-cluster 
of C born in [0,t]. This is done in the following way. First, those C G A[o,t](C) that have 
no ancestors are kept. Then we look to the remaining cylinders in A[o,f](C) and erase those 
that have a kept cylinder in its first generation. We repeat these two steps for the cylinders in 
A[o,t](C) that have not already declared to be kept or erased, and continue in this way until 
we reach C. The end result is a partition of A[o,j](C) in two subsets formed, respectively, by 
kept and erased cylinders. In particular, the subset to which C belongs decides its status. In 
fact, under (|2.7|) one can prove that all cylinders have a finite number of ancestors born in the 
interval [0,t], and, thus, the process r/j can be constructed following, BO-cluster by BO-cluster, 
the steps of the finite case ( ^.1[ ). In the next section, we sketch the proof of the finiteness of 
the number of ancestors and give further details of the construction. 

It is natural, and it turns out to be convenient, to extend the notion of ancestors of a cylinder 
to that of ancestors of a space-time point (x,t), x G Z'^ and t G M: Let the first generation 
of ancestors of (x, t) be the set of cylinders in C whose basis contains x and are alive at time 
t. The nth generation of ancestors of (x, t) is then formed by the (n — l)-th generation of 
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ancestors of the cylinders in the first generation. The union of all the generations of ancestors 
is the BO-cluster A{x,t) of {x,t). More generally, the set of ancestors of T C Z'^ at time t is 
defined by 

A(T,t)= [jAix,t). (3.2) 

4 Existence of /i and exponential convergence 
4.1 Backwards percolation 

To perform the construction described in the previous section, every cylinder C G C must have 
a finite number of ancestors. If this is the case, we say that there is no (backwards oriented) 
percolation in C. Hence, if with probability one there is no backwards percolation, the double 
infinite construction holds and we have a process {rit)teR that is time-translation invariant. The 
marginal distribution of rjt does not depend on t and it is called fi. By construction fi is an 
invariant measure for rjt. This shows the existence of /i in a constructive way. In contrast, the 
existence of fi under (|2.7|) uses a fixed point theorem. 

The condition jS > 13m implies that there is no percolation with probability one. This is 
shown by dominating the number of plaquettes in the bases of the cylinders in a BO-cluster by 
a branching process. The number 

{d-l)\p (4.1) 

is an upperbound on the mean number of branches of the process. That is, the mean number 
of plaquettes born from the branching (= incompatible contours) of each single plaquette. The 
process is subcritical if this number is less than one, thus the condition j3 > (3m- This argument, 
sketched in Section | below is inspired by Hall (1985), who dominated a continuum percolation 
process by a branching process. We sketch this domination in Section ^ below. If there is no 
percolation, the number | A(x, t) \ is finite for all (x, t). As a consequence, there exists a function 
$ : (/, A(Supp {f),t)) $(/, A(Supp (/), t) such that for any / with finite support 

/(r^,)='^'(/,A(Supp(/),t)). (4.2) 
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For instance, to decide wiiether a contour 7 is present at time t it suffices to look at the BO- 
cluster of {x, t) for some a; G 7. The function $ is the one that decides which cyhnders are kept 
and indicates the presence/ absence of 7 at time t. 

4.2 Time length and space width of the BO-cluster 

Most of the stated properties — uniqueness of /i, exponential clustering and finite-volume effects, 
and exponential convergence to equilibrium of the loss network — follow from the observation 
that for P > Pm both the time length and the space width of the BO-cluster of any given site 
decay exponentially. More precisely, let us introduce Proj (A(a;,t)) C Z*^, the spatial projection 
of the BO-cluster, defined as the set of points in IJ^ belonging to the basis of some cylinder in 
the BO-cluster: 

Proj(A(x,i))= U {^£7} (4.3) 

7ea(a;,t) 

where a(x, = {7 : (7, Tjt(7), Sk{'y)) £ A(x, t) for some k}, is the set of bases of the cylinders of 
the BO-cluster of (x, t). The cardinality of this set will be bounded by the cumulative number 
of points: 

||A(a:,t)||= Yl 1^1- (4.4) 

Indeed, it is clear that 

|Proj(A(a:,t))|<||A(a:,t)||. (4.5) 

We then have: 

1. Let E2{t) be the set of C for which the BO-cluster of (0, 0) has time- length larger than t: 
E2{t) = {C : C is alive at time - t for some C e A(0, 0)} . (4.6) 
Then, for (5 > (5m 

P(£;2(i)) < Ml e-*(i-(''-^)^/5) (4.7) 

with Ml > 0. 
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2. Let E^i^w) be the set of C for which the projection of the BO-cluster of (0,0) is not 
contained in [—w,wY: 



E,{w) = |C : Proj (A(0,0)) ^ [-i/;,^'}- 



Then, for /? > /?. 



HEsiw)) < M2e 



-M3W 



where M2, M3 are as in ( 2.17] ). 



The proof of (|4.7| ) and (|4.9|) are sketched in Section ||. To prove 
a continuous-time branching process. On the other hand, to prove 
by the total population of a branching process. 



(4. 



(4.9) 



we dominate A(0, 0) by 
we dominate ||A(0,0)|| 



4.3 Proof of Rl and R2 

The exponentially fast time-convergence ( |2.16| ) is a consequence of ( ^^ ) and (|4.9|). We use 
the same Poisson marks to construct simultaneously the stationary process rjt and a process 
starting at time zero with an arbitrary initial configuration ^. The second process is called 
C,t, where ^0 = ^- The process C,t ignores the cylinders in C with birth-times less than and 
considers cylinders in C(0) with basis in ^ and birth-time zero. The process rjt ignores the 
cylinders in C(0). Hence for any 7^0, 

\Vt{l)-Ul)\ < l{A(0,t) ^ A(0,t)} (4.10) 

where A(0,t) is the cluster constructed in C[o,t] U C(0). In FFG it is shown, via a coupling 
argument, how ( ]4.8| ) and ( |4.9| ) imply that the expectation of the right hand side of (|4.10| ) decays 
as exp(— Mot). The exponential decay of length and width of the cluster implies exponential 
decay of the probability that it contains a cylinder of the initial configuration C(0). 
The uniqueness of fi follows immediately from ( |2.16| ). 

Reversibility follows from the facts that (rj^) converges in distribution to (rjt), fi^ converges 
to fi and fi^ is reversible for rj^. From the construction, under /3 > Pm, it is possible to show 
that {r]^) converges almost surely to (r/t). Some details are given in the next sections. 
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5 Space-time mixing and the central limit theorem 
5.1 The key facts 

The mixing properties of the measure n are a consequence of the following space-time mixing 
properties of C. 

• Let / be a function depending on contours contained in a finite set A. Let rjf^ be the loss 
network process constructed in A. Then 

|E(/M)-E/(r/o^)| < 2||/|Up(A(Supp(/),0)^A^(Supp(/),0)). (5.1) 

where A'^(Supp is the cluster of (Supp constructed with cylinders in 

= {(7, n(7), Skil)) G C : 7 C A, A; G Z}, (5.2) 

the subset of cylinders whose basis is in A. 

• For arbitrary measurable functions / and g, 

mfiVo)9ivo)) - ^fiVo)^9iVo)\ 
< 2 ||/||oo||^||ooP(c" n C'V for some C G A(Supp (/), 0) and C" G A(Supp (g), 0) 

(5.3) 

where A(Supp ((?), t) has the same distribution as A(Supp (^f), t) but is independent of 
A(Supp(/),t). 

The proof of (|5 . 1| ) follows rather straightforwardly from the space-time construction. Using 
(13) we get 

f{Vo)-f{Vo) = [$(/,A(Supp(/),0))-$(/,A^(Supp(/),0)) 

X l{A(Supp (/), 0) ^ A^(Supp (/), 0)} . (5.4) 

As, by definition, |$(/, A(Supp (/), t))! < ||/||oo, taking expectations and absolute values in 
(U we get (|51|). 

The proof of ( |5.3| ) is similar in spirit but requires a somewhat more delicate argument based 
on the coupling of two continuous-time versions of the backwards percolation process. See 
details in FFG. 
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5.2 Proof of R3 and R4 

To prove the finite-volume effects ( |2.17| ) we use the space-time representation ( |3.1| ) and get 

/i/-/iV = E/(r7o)-E/(r/o^). (5.5) 

By ( |5.1| ) it is enough to bound 

p(a(Supp (/), 0) + A^(Supp (/), 0)) , (5.6) 
which as in ( |4.9| ) is bounded by 

M2 e-*^^'^^^'^^) . (5.7) 

xGSupp (/) 

This proves the decay stated in (|2.17|) . 

The proof of exponential mixing ( p. 181) is similar but using instead the bound (|5 



While we have not yet done a careful study, we believe that (|2.17|) and (|2.18|) lead to 



sharper inequalities than those obtained via the use of "duplicated variables" [von Dreifus, 
Klein and Perez (1995), Bricmont and Kupiainen (1996)]. The reason is that clusters formed 
by superposition of two systems of contours have larger probabilities of intersection than our 
single-system clusters. 

5.3 Proof of the central limit theorem 

We use the results for stationary mixing random fields of Bolthausen (1982). Let = r^f ■ 
By hypothesis, ||X^||2+<5 < oo. Under this conditions, Bolthausen (1982) shows that if 

oo 

Y,n'-\a,,^{n)Y'^'^'^ <oo (5.8) 

n=l 



then D < oo and (|2.19|) holds. Here 02,00 (^) measures the dependence between functions 
depending on the sigma algebra generated by Xq and Xy and the sigma algebra generated by 
{Xx : x G A} for |A| = 00 and min{|a;|, \y — x\ : x E X} > n. In FFG we use (|2.18|) to show 



that 



«2,oo(n) < (M2)2|Supp(/)| (5.9) 

||/|>n-2|Supp(/)| 
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Hence, a2,oo{n) decreases exponentially fast with n. This shows the central limit theorem. 



6 Length and width of the BO-cluster 

To conclude, let us sketch the arguments behind the bounds ( ^?7| ) and ([4.91). In both cases we 
rely on dominating branching processes. 

6.1 Time length 

To show ([4.7|) we consider a continuous time multitype Markov branching process ht on N^. In 
this process, each contour 7 lives an mean-one exponential time after which it dies and gives 
birth to kg contours 6, 6 E Q, with probability 

n-*'"'r'7-^''^° (6.1) 

for kg > 0. These are independent Poisson distributions of mean /i(7, 6) = 1(7 fl 7^ 0} e^^'^L 
Fix bo(7) = \{k : (7, (7), 5*^(7)) is alive at time 0}| 1(7 3 0} and zero otherwise. Under 
this initial condition it is possible to couple (b()(>o and A(0, 0) in such a way that 

E,{t) C \ j2^,{e) = 0]. (6.2) 

e 

Using the backwards Kolmogorov equation for Rf = E^gbt(^), one can show that 

> 0) < e""^-!)^^-!)*. (6.3) 

9 

6.2 Space width 

Define a Galton- Watson branching process Z„ G N as follows. Let Fj" be i.i.d. non negative 
integer valued random variables with the same distribution as 

F:=5^(rf-l)|7|X, (6.4) 

790 
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where are independent integer valued random variables with Poisson distribution of mean 
e-^lTl. Define Zq = I and 

^".+i = E^" (6.5) 

i=l 

(with the convention Y^i=i = 0)- is possible to couple the BO-cluster A(a;, t) and {Zn)n>o 
in such a way that the number of plaquettes in the bases of the cylinders in the nth generation 
of ancestors of (x, t) is less than or equal to Zn- 



||A(x,t)|| < (6.6) 

n>0 

Hence, to show ( [4.9|) it suffices to prove 

P(Z >k)< M2e-^'^ (6.7) 

where Z = ^„>o-^n. Call F{b) the generating function of Z, we will prove that if /3 > Pm, 
b = sup{6 : F{b) < oo} > 1. 

The generating function of Y is given by 

/(a) = Ea^ = JjEa^'^-i)!^!^- = exp (j^e-^^^^ia'^'''^^^'^ - 1) ] . (6.8) 

790 \790 / 



The radius of convergence of /(a) is given by exp(/5 — /5p), where (3p is defined in (|2.8| ). For 
> (3m{> l3p), the radius of convergence is strictly larger than 1. The mean number of 
offsprings KZi is given by 

f{a)l^, =EY = {d-l)Y, ble-'^l^l < 1 (6.9) 

790 

for j3 > (3m- Hence, our branching process is subcritical, 

1 = /(I) and X = f{x) implies a; > 1 (6.10) 

i.e. the smallest solution of the equation x = /(x) is 1. 

By (13.3) of Harris (1963) F{b), the generating function of Z, must satisfy the equation 

F{b) = bf{F{b)). (6.11) 
20 



The largest solution of this is 

h = -a/f{a) (6.12) 



where a is the solution of 



a 

In this case, it is easy to see that 



/'(a) = (6.13) 



/'(a) = _ 1) ^ |7|e-^l^la(^-i)l^l (g^^^^ 

(X 

730 



and a is the solution of 

^ |7|e-^l'^la(''-^)l'^l = l/{d - 1) (6.15) 

790 

which gives us 

a = e('^-^M)/(d-i) (Q_1Q) 

Therefore, 

I = expj^l^ + Y^e^Mb\/(d^i)^i _ e{/5-/3M)l7l/(<i-i))|. (6.17) 

^ ' 730 

By exponential Chebichev, fixing M2 = IE6^ and M3 = log 6, we get ( |6.7| ). 
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